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1.  Introduction 


In  this  paper  we  review  some  results  concerning  the  total  temperature  sensitivity  of 
constant-temperature  anemometers  in  supersonic  flows.  In  particular,  the  non-linear  response  at 
low  overheat  ratio,  first  presented  by  Smits  &  Perry  [1]  and  Smits  et  al.  [2]  is  discussed. 

2.  CTA  sensitivity  to  total  temperature  fluctuations 


At  supersonic  velocities,  a  CTA  is  sensitive  to  both  mass  flow  (pu)  and  total  temperature  (To) 
fluctuations  [3].  Therefore,  a  first  order  small  perturbation  analysis  of  the  output  voltage  e  yields: 

e'=  kpu-(puy+kT()-T0'  ,  (1) 

where  kpu  and  kTo  are  the  anemometer’s  sensitivity  coefficients  to  mass  flow  and  total 
temperature  fluctuations,  and  the  prime  denotes  fluctuating  quantities.  When  the  overheat 
parameter  r  =  (Tw  -  ijT0)  /  T0  and  the  total  temperature  are  constant  (q  is  the  wire  recovery 
factor),  kpn  is  obtained  by  plotting  the  anemometer  output  voltage  as  a  function  of  the  mass 
flow  following  the  semi-empirical  relation  [2]: 


which  yields: 


e2  =L(T0,T)+M(T0,T)-(pu)n, 
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where  Sp,,  is  the  non-dimensional  mass  flow  sensitivity.  It  can  be  seen  in  (3)  that  kpu  and  Sp,, 

are  function  of  the  three  variables  i,  T0,  and  pu.  However,  L  is  usually  small  so  that  S/M  can  be 
approximated  by  Spi,~n/ 2. 

Since  an  independent  variation  of  the  total  temperature  with  constant  mass  flow  is  often 
difficult,  kTo  is  usually  deduced  from  the  experimental  value  of  kpu  obtained  using  the 

precedent  calibration  procedure.  Indeed,  differentiating  (2)  with  respect  to  T0  leads  to  the 
following  relation  for  the  total  temperature  sensitivity  [2,  4]: 
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where  a,  b,f(r),  and  g(r)  are  defined  in  [2].  It  can  be  seen  that  STo  doesn’t  depend  directly  on  Tn 
and  pu.  These  dependences  are  hidden  in  Spil  and  t.  In  particular,  when  Spu  is  assumed  constant, 
STo  depends  only  on  t. 

This  definition  of  kj  takes  into  account  the  variation  of  the  flow  properties  with  total 
temperature  (constants  a  and  b)  and  the  variation  of  heat  transfer  with  overheat  ratio  (functions/ 
and  g ).  When  these  parameters  are  neglected  and  when  the  recovery  factor  q  is  assumed  to  be 
unity,  equation  (4)  degenerates  in  the  following  simple  relation,  which  is  generally  used  in  sub¬ 
sonic  flows  [5]: 
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Figure  1  shows  a  plot  of  STo  for  the  two  cases  (subsonic  and  supersonic),  using  the 
calibration  data  of  Kovasznay  [3]:  77=0.95,  a=b= 0.768,  /(r)=g'(r)=l-0.18r,  n=  0.5, 

5^=77/2=0.25.  The  diagram  in  the  inset  show  the  same  curves  with  a  zoom  on  the  y-coordinate: 
it  shows  that  equation  (4)  leads  to  negative  values  of  STo  when  x  is  high  enough  (t  >  1).  This  is 
physically  impossible  and  the  whole  relation  breaks  down  at  high  overheat,  because  of  wire 
burnout.  The  problem  arising  at  low  overheat  ratio  is  clearly  seen  in  figure  1 :  a  small  variation 
of  T0  produces  a  variation  of  x  if  Tw  is  maintained  constant  ;  therefore,  since  STo  is  a  strong 
function  of  x  at  low  overheat,  the  instantaneous  sensitivity  can  differ  significantly  from  its  value 
at  the  average  total  temperature,  thus  compromising  the  linear  based  modal  analysis  [2]. 

3.  Non-linearity  of  the  response  to  total  temperature  fluctuations 

A  more  quantitative  approach  to  this  issue  can  be  obtained  by  developing  the  anemometer’s 
response  to  total  temperature  fluctuations.  If  the  mass  flow  over  the  wire  is  supposed  constant,  a 
variation  AT0  of  the  total  temperature  results  in  a  variation  Ae  of  the  output  voltage: 

Ae  =  kTQAT0+R(AT0)  ,  (6) 

where  R  consists  of  terms  of  order  greater  than  one.  If  the  system  were  perfectly  linear,  k  T 

would  be  constant  and  R  would  be  zero.  This  is  in  reality  not  the  case,  but  when  A Tn  is  suffi¬ 
ciently  small,  R  is  negligible  and  a  local  linearization  approach  (ie  a  first  order  analysis)  may  be 
sufficient.  Flowever,  since  k T  is  (like  STo)  a  strong  function  of  x  (and  consequently  of  T0)  at  low 

overheat,  higher  order  terms  may  have  an  influence  even  for  reasonably  small  variations  of  T0, 
thus  compromising  the  first  order  analysis.  Equation  (6)  can  be  written  in  the  form 
Ae  =  kjQ  Ar0(l  +  P) ,  where  P  =  R(AT0)  /  kT  AT0  is  a  measure  of  the  error  induced  by  neglect¬ 
ing  a  finite  value  of  R.  When  P  is  positive,  the  total  temperature  variation  is  overestimated. 

If  a  known  heat  transfer  relation  is  assumed,  P  can  be  obtained  for  a  given  value  of  A T0  by 
computing  the  actual  variation  Ae  of  the  anemometer  output  voltage  using  equation  (2),  as  well 
as  the  temperature  sensitivity  k T  using  (4)  (algebraic  relations  for  L  and  M  are  given  in  refer¬ 
ence  [6]).  Results  showing  the  variation  of  P  as  a  function  of  the  overheat  ratio  are  presented  in 
figure  2.  As  before,  the  calibration  data  of  Kovasznay  has  been  used.  The  curves  were  obtained 
by  computing  the  difference  Ae  =  e(T0  +  AT0  /  2)  -  e(T0  -  ATq  /  2)  and  comparing  it  with 
kT  (To)-  ATq  to  obtain  R  and  then  P.  The  mean  total  temperature  Tq  is  293K  and  two  values  of 

A  To  have  been  chosen:  AT0=48K  corresponds  to  a  variation  of  +/-  2  a  where  <r=4%  is  a  typical 
total  temperature  RMS  fluctuation  in  a  supersonic  flat  plate  turbulent  boundary  layer  [2].  In 
this  case,  P  stays  lower  than  1%  only  in  a  limited  domain  of  overheat  ratios  comprised  between 
7=0.3  and  7=0.9.  A  full  modal  analysis  down  to  t<0.1  appears  therefore  to  be  impossible,  as 
suggested  by  Smits  &  Dussauge  [7].  On  the  other  hand,  for  A7’0=12  K  (which  corresponds  to  a 
variation  of  +/-  2a  where  cr=l%),  the  error  appears  to  be  acceptable  down  to  7=0.05. 

According  to  this  analysis,  investigation  of  turbulent  flows  involving  very  low  total  tem¬ 
perature  fluctuations  is  therefore  possible  with  a  constant  temperature  anemometer.  This  is  for 
example  the  case  in  the  free-stream  of  supersonic  wind  tunnels,  where  very  small  reminiscent 
“entropy  type”  disturbances  may  be  superimposed  to  acoustic  waves  (the  same  conclusion  has 
been  reached  in  reference  [8],  albeit  with  a  slightly  different  method).  These  results  show  how¬ 
ever  that  when  the  total  temperature  fluctuation  reaches  several  percent,  a  small  perturbation 
analysis  of  the  output  voltage  fluctuation  is  insufficient,  as  argued  by  Smits  &  Dussauge  [7]. 
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4.  Influence  of  the  feedback  system 

The  analysis  presented  in  sections  2  and  3  were  made  assuming  that  the  feedback  system 
imposes  a  constant  value  of  the  wire  temperature  Tw.  In  reality,  the  wire  resistance  has  to  vary 
by  a  small  amount  to  give  the  amplifier  a  finite  feedback  signal,  so  that  the  assumption  of 
constant  wire  temperature  is  only  an  approximation. 

The  static  analysis  of  a  non-ideal  CTA  has  been  performed  by  Perry  &  Morrison  [9]:  in  essence, 
the  operating  resistance  Rw  of  the  wire  is  given  by  the  intersection  of  two  equations: 

-  The  equation  relating  the  output  voltage  to  the  flow  parameters  ;  this  is  the  square 
root  of  equation  (2): 

e  -  tJl(Rw)  +  M(RW)-  (pu)n  (7) 

-  The  equation  obtained  by  an  analysis  of  the  electrical  circuit  [9]: 

(Ra+Rw)-(Rb+Rc)-Eqi 

e  = -  (8) 

(Ra  +  Rw)'  (Rb  +  Rc)  +  Kdc  '(Rw  '  Rc  ~  Ra  '  Rb) 

where  Ra,  R/„  and  Rc  are  the  resistances  of  the  Wheatstone  bridge,  Kdc  is  the  amplifier’s  DC 
gain,  and  Eqi  is  the  anemometer’s  offset  voltage1. 

The  wire  operating  point  is  best  represented  graphically  by  the  intersection  of  the  curves 
(7)  and  (8)  on  a  ( e ,  Rn)  plane:  this  is  presented  in  figure  3  for  an  overheat  parameter  of 

approximately  t~0.3  (the  feedback  amplification  is  7f/,=5000).  It  can  be  seen  that  the  wire 

resistance  is  always  higher  than  the  resistance  corresponding  to  a  perfectly  balanced  bridge,  and 
that  it  depends  on  the  offset  voltage  Eqi\  the  higher  the  offset  voltage,  the  higher  the  wire 
resistance.  In  particular,  it  is  important  to  notice  than  a  very  high  DC  gain  is  not  sufficient  to 
balance  the  bridge.  As  shown  by  Perry  &  Morisson  [9],  the  Wheatstone  bridge  would  be 
perfectly  balanced  only  for  Kdc  — >  oo  and  Eql=  0.  For  the  present  case,  taking  values  of  Kdc  larger 
than  5000  has  a  negligible  influence  on  the  operating  point.  The  offset  voltage  has  a  large 
influence  on  the  anemometer’s  dynamic  performance  and  typical  values  of  Eqi  range  between 
0.1  mV  and  30  mV,  depending  on  the  anemometer’s  characteristics. 

Since  the  curve  corresponding  to  equation  8  is  not  vertical,  the  wire  resistance  varies 
slightly  with  the  temperature  and  the  actual  sensitivity  of  the  real  anemometer  is  therefore  lower 
than  the  one  obtained  by  considering  a  sensor  maintained  at  a  constant  resistance.  The  deviation 
from  an  ideal  anemometer  is  a  function  of  the  offset  voltage  and  of  the  overheat  ratio.  At  very 
low  overheat,  the  temperature  sensitivity  decreases  drastically,  as  shown  in  figure  3  as  well  as 
in  the  figure  3  of  reference  [1]. 

Consequently,  the  total  temperature  fluctuation  measured  with  a  real  system  by  using  the 
sensitivity  obtained  with  equation  4  or  5  will  be  under-estimated,  particularly  at  low  overheat 
and  high  offset  voltage.  This  is  illustrated  in  figure  4,  where  the  error  P  is  computed  by  taking 
this  effect  into  account.  For  the  seek  of  simplicity,  fluid  property  variations  are  ignored  (i.e., 
a=b= 0  and ./(r)=g(r)=0). 

in  opposition  to  the  results  obtained  in  section  3,  P  appears  to  be  independent  on  the  fluctua¬ 
tion  level  as  long  as  Eqi  is  not  too  small.  For  Eqi= 0,  the  system  behaves  like  a  perfect  constant  tem¬ 
perature  bridge  and  P  has  the  same  value  as  in  figure  2.  For  Eqi> 0,  the  total  temperature  fluctua¬ 
tions  appear  to  be  underestimated,  and  the  relative  error  is  significant  even  at  high  overheat  ratio. 

It  should  be  noted  that  the  same  phenomenon  happens  for  the  sensitivity  to  mass  flow 
fluctuations.  The  real  anemometer  sensitivity  is  lower  than  the  ideal  one  because  of  the  finite 
offset  voltage.  In  this  case,  however,  the  deviation  is  taken  into  account  during  calibration  since 
the  mass  flow  sensitivity  is  obtained  by  recording  the  output  voltage  of  the  real  anemometer  as 
a  function  of  the  mean  mass  flow,  so  that  the  operating  point  follows  curve  (8)  during  the 


1  A  similar  expression  for  a  CTA  with  offset  current  has  been  derived  in  reference  [8]. 
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Fig.  1:  CTA  Total  temperature  sensitivity  Fig.  2:  Response  non-linearity  (ideal  CTA) 

(a)  supersonic  case  (eq.  4) 

(b)  subsonic  case  (eq.  5) 
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Fig.  5:  Error  in  the  overheat  estimation 

calibration  procedure.  As  pointed  out  by  Perry  &  Morrison  [9],  the  effect  of  wire  resistance 
variation  with  the  mean  mass  flow  is  therefore  to  change  the  coefficient  M  in  equation  (3)  so 
that  the  gradient  of  the  calibration  curve  deviates  from  the  ideal  constant  temperature  value.  The 
measured  slope  M(T0,  t)  is  smaller  than  the  one  which  would  be  obtained  with  an  ideal 
anemometer  and  account  for  the  decrease  in  kpu,  as  can  be  seen  by  considering  equation  (3). 
This  has  however  no  effect  on  the  total  temperature  sensitivity,  as  shows  equation  (5). 

These  results  show  that  a  derivation  of  kT  assuming  a  constant  wire  temperature  can  lead  to 
serious  errors  when  total  temperature  fluctuations  are  measured,  regardless  of  the  level  of  these 
fluctuations.  Strictly  speaking,  these  errors  are  not  due  to  the  non-linear  heat  transfer  relation,  but 
to  the  inability  of  the  feedback  loop  to  maintain  a  constant  wire  temperature.  To  avoid  these  er¬ 
rors,  a  direct  calibration  of  the  CTA  at  different  total  temperatures  should  be  performed.  This  re¬ 
quires  a  facility  allowing  independent  variation  of  the  mass  flux  and  the  total  temperature.  An¬ 
other  alternative  would  be  to  use  an  anemometer  with  frequency  dependent  gain,  where  the  offset 
voltage  has  no  direct  influence  on  the  dynamic  behaviour  and  can  be  maintained  very  small  [10]. 
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5.  Influence  of  the  offset  voltage  on  the  overheat  ratio 

It  has  been  showed  in  the  last  section  that  the  wire  resistance  in  a  real  CTA  is  always 
slightly  higher  than  the  resistance  corresponding  to  a  perfectly  balanced  bridge: 
(Rw)\n\m<x=(Ra'Rb)/Rc-  This  means  that  the  exact  wire  overheat  depends  on  the  flow  conditions 
and  Eqi  and  is  not  known  a  priori.  Very  often,  the  exact  determination  of  the  wire  overheat  is 
not  of  principal  importance  if  the  measurements  are  performed  at  the  same  settings  as  during 
calibration,  but  when  the  fluctuation  diagram  technique  is  used,  the  wire  overheat  becomes  an 
important  parameter  and  its  determination  should  be  made  accurately.  This  is  especially  the 
case  when  the  anemometer  allows  an  offset  variation  for  frequency  optimisation  purpose:  since 
measurements  at  supersonic  velocities  require  a  careful  adjustment  of  the  dynamic  response,  a 
change  of  Eqi  can  usually  not  be  avoided  during  the  tests. 

The  relative  error  which  is  made  by  assuming  a  perfectly  balanced  bridge  is  presented  in 
figure  5.  The  real  (or  “effective”)  overheat  parameter  is  denoted  Teff  whereas  the  overheat 
corresponding  to  a  perfectly  balanced  bridge  is  written  xbal.  Since  kT  is  a  strong  function  of  the 
overheat,  neglecting  this  error  can  lead  to  an  overestimation  of  the  sensitivity.  To  overcome  this 
problem,  the  authors  propose  to  record  the  anemometer  output  voltage  and  the  wire  voltage 
simultaneously,  in  order  to  be  able  to  compute  the  effective  wire  overheat. 

6.  Conclusion 

It  has  been  showed  that  the  non-linear  relationship  between  total  temperature  variations 
and  the  output  voltage  of  a  CTA  can  induces  large  errors  at  low  overheat  ratios,  as  already  dis¬ 
cussed  by  Smits  &  Perry  [1]  and  Smits  et  al.  [2].  This  is  however  restricted  to  relative  fluctua¬ 
tions  of  several  percents,  typical  of  fully  developed  turbulent  flows.  For  a  lower  fluctuation 
level,  like  in  the  free  stream  of  supersonic  wind  tunnels,  a  first  order  approximation  seems  accu¬ 
rate  enough  to  perform  a  full  modal  analysis.  Of  more  concern  is  the  overestimation  of  total 
temperature  sensitivity  that  occur  when  performing  only  a  direct  mass  flow  calibration  and  de¬ 
riving  a  relation  assuming  constant  wire-resistance,  since  it  appears  to  be  independent  of  the 
fluctuation  level. 
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